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This paper analyzes different methods of adjusting the sampling time 
for detecting synchronous binary data, based on properties oj the random 
data signal itself. The static error and the variance of the jitter of the 
resultant sampling instant are calculated where the effects of frequency 
offset, additive noise, signal overlap, and jitter of the reference source 
are included. 

The threshold crossing timing recovery system adjusts the sampling time 
in response to the times at which the data signal crosses the amplitude 
threshold. The sampled-derivative system tises the time derivative of the 
signal at the sampling time to adjust sampling phase. It is shown that 
both systems lead to approximately the same amount of jitter in the presence 
of noise and signal overlap for a given bandwidth of the control loop. 

An improved timing recovery system is presented which is constructed 
by adding correction signals to the sampled-derivative system. This system 
accounts for intersymbol interference in a manner that tends to set the 
sampling time at the point of maximum eye opening, where the error 
probability is minimum for the most adverse message sequence. 

I. INTRODUCTION AND SUMMARY OF RESULTS 

In synchronous polar binary data transmission, information is sent 
by serially transmitting either a basic signaling waveform or its negative 
at fixed time intervals. Modulation may be used to better fit the signal 
to the channel. At the receiving end, the signal is demodulated and 
filtered. The resultant baseband signal is sampled periodically, and 
the polarities at the sampling instants determine the output data. 
The choice of sampling time is critical for minimizing the error prob- 
ability due to intersymbol interference and noise, particularly when 
the signal has been subjected to sharp cutoff filtering. The sampling 
time is best set by using some properties of the data signal itself. 
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The problem of timing is particularly acute in pulse code modulation 
(PCM) systems, where the accumulation of jitter in a long chain of 
regenerators frequently limits the allowable; length of such systems. 
For this reason, previous studies of timing recovery have concentrated 
on PCM applications. 1 " 4 The use of a tuned circuit as the memory 
element is generally assumed, since this is commonly employed in 
PCM repeaters. 

This paper will concentrate on timing recovery for data transmission 
applications. The effects of multiple regeneration will not be considered. 
The recovery of timing will be accomplished by a feedback control 
system, such as a phase-locked loop. Different methods of generating 
the error signal for the control loop will be compared. 

The received signal, after demodulation and filtering, is of the form 

s(f) = E a k f(t - 0T - kT) + n(t) (1) 

where \a k \ is a set of independent random variables, each equal to 
+ 1 or —1 with equal probability. This may be assured by the use 
of a scrambler if the data source is itself not random. The basic signal- 
ing waveform is j(t). The abscissa of f(t) will be adjusted for each system 
to be studied so that the desired sampling time of f(t) is t = 0. The 
quantity is an unknown fractional time delay. Since we are not 
concerned with absolute time delay between transmitter and receiver, 
we will assume | /3 1 ^ %. The additive noise is n(t). 

The sampling wave which determines the times at which s(t) is 
sampled may be represented by 

q(t) = f) 8(i 1-tiT ~yT), (2) 

n = — oo 

where 7 is a phase that is generally time varying. 
The output data is determined by 

d„ = sgn s(7iT + yT), (3) 

where sgn v = v/\v\. Then 

a n = sgn [a n j(yT - &T) + £ a n - k j(kT + yT - (3T) + n(f)]. (4) 

For simplicity, the argument of the noise term is not made explicit 
since it is of no consequence. Assuming that 1(yT — &T) is positive, 
then d„ will agree with a n provided that 

-ajE a n - k j(kT + yT - 0T) + n(f)] < f(yT - 0T). (0) 
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It is readily seen that the probability that (5) holds depends strongly 
on 7 — /S. For each timing recovery system, f(t) is denned so that the 
desired value of y — /3 is zero. 

The principal part of a timing recovery system is a phase detector 
which examines s(t) and q(t) and attempts to generate an error signal 
proportional to /3 — 7. It is not possible to determine /3 exactly from 
s(t) since the signaling waveform and the noise are unknown. This 
paper will consider different methods of forming this estimate of /3 — 7. 

Another essential component of the system is the reference source 
which is used to generate the sampling wave. Its phase or frequency is ad- 
justed by the error signal in order to form a sampling wave of the proper 
phase. The error signal may be filtered prior to its use in shifting the 
phase of the reference source. The reference source may be a local 
oscillator whose natural frequency is set as close as possible to the 
bit rate. The reference source may instead be derived from transmitted 
pilot tones, in which case its frequency is exact, but it might have 
phase jitter of its own due to channel noise. 

Section II describes and analyzes a timing recovery system which 
uses a threshold crossing phase detector. This detector generates an 
error signal each time the signal crosses zero. The amplitude of this 
signal is proportional to the difference between the time of occurrence 
of the zero crossing and the time of the nearest sampling pulse, dis- 
placed by half a bit period. This system tends to choose a sampling 
instant which is midway between the mean transition times. 

The sampled-derivative phase detector is discussed in Section III. 
This device generates an error signal during each bit interval which 
is proportional to the time derivative of the signal at the sampling time 
multiplied by the signal polarity at that time. The sampled-derivative 
timing recovery system attempts to set the sampling time to coincide 
with the peak of j(J). 

The analysis shows that the performance of these two systems is 
very similar for a given open loop gain function of the control system, 
G(o)). Approximations are made based on the assumption that the 
phase error is small and that G(u) is a narrowband low-pass function 
compared with the bit rate. 

The systems fail if (7(0) is finite and the reference source does not 
agree exactly in natural frequency with the bit rate. If the reference 
source has the correct frequency and a phase 0, then a static phase error 

6 = T+W) (6) 

results in the sampling wave. 
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Much better performance can be achieved if G(u>) has a pole at the 
origin. In this case the system is insensitive to the phase of the reference 
source. In the presence of a frequency error, A/, the static error is 

_1 + G(oj)J u _o 

In addition to any static error, the sampling time will also jitter 
about its mean value. The variance of this phase jitter is the sum of 
several components, each due to a different cause. Jitter is produced 
by jitter in the reference source, by additive noise and by signal overlap. 
In the case of the threshold crossing system, jitter is also introduced 
whenever there is a static error. 

If the reference source has a jitter whose power spectral density is 
S,(w), then the output fractional jitter will have a variance equal to 

°"s = -] n |i + G»| 2d£0 - (8) 

This indicates that high-frequency noise components must be removed 
from the reference source prior to its use for timing recovery. 
The jitter produced by the additive noise is 

2 _ 2[/? n (0) - R n (T)\ m 

<r N - r[/ '(_ r/2 ) -f(r/2)f B| w 

for the threshold crossing system. /?„(/) is the autocorrelation of the 
noise and co, is the noise bandwidth of the closed control loop. 



_ i r G(u) 

0,1 " X i„ I 1 + C(co) 



du. (10) 



For the sampled-derivative system, the noise leads to jitter variance 



2 ^ R'„'{0) 



(11) 



In typical data transmission systems, (9) and (11) are similar in 
magnitude, and not very sensitive to the shape of /(/) if the noise is 
similarly filtered. 

The jitter variance due to signal overlap is of the form 

4 = AvJ + AMT)\ (12) 

where 

3 i r . I gfe) 

W2 = r h u I 1 + G(o) 



do>. (13) 
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If b) x T and <a 2 T are comparable and much less than unity, then the 
first term is usually much larger than the second. 
For the threshold crossing system, 

•[2/(fcr + T/2) + K-kT + T/2) - K-kT - T/2)] (14) 
and 

Az = 2T\]\-T/2) - f(T/2)) 2 \-*fW 2 ) 
+ 22 f(kT + T/2)f(kT - T/2) - f) k 2 \j(kT + T/2) 

k-— oo *•=•— oo 

- j(kT - T/2)][f(-kT + T/2) - j(-kT - T/2))}- (15) 

For the sampled-derivative system, 

Ai = tTW) t ?L www + n-kT)] (i6) 



and 



A 2 = 2rj ,} 2{0) JC k 2 r(kT)f(-kT). (17) 

In both cases, A, = if /(0 is an even function, so the timing recovery 
systems are very sensitive to asymmetry of the basic signaling wave- 
form. The jitter variances are again comparable for both systems. 
As may be expected, the jitter increases considerably as the filter used 
to shape f(t) is made sharper. 

There is an additional jitter component for the threshold crossing 
system whenever there is a static error. Its variance is given by 

a] = e 2 ^T. (18) 

An example is provided in Section IV. A typical data transmission 
system using a distorted signal is studied so as to illustrate the mag- 
nitudes of the above quantities and to indicate the narrowness of loop 
bandwidth required for satisfactory performance. 

In this example it is also seen that neither the threshold crossing 
timing recovery system nor the sampled-derivative system chooses a 
mean sampling time which is very near to the time at which the eye 
pattern has its maximum opening. 
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Section V describes an improved timing recovery system whose 
mean sampling time coincides with that of the maximum eye opening. 
This system is constructed by adding correction signals to the sampled- 
derivative system in order to account for the effects of intersymbol 
interference on the mean sampling time. 

Finally, an outline of some extensions and modifications of these 
timing recovery systems is presented in Section VI. 



II. THE THRESHOLD CROSSING SYSTEM 

Most timing recovery systems make use of the instants at which 
the data signal crosses the threshold to alter the phase of the sampling 
wave. A block diagram of a typical threshold crossing timing recovery 
system is shown in Fig. 1 . 

The principal part in this system is the threshold crossing detector. 
This device generates an error pulse each time the signal crosses zero. 
The amplitude of the error pulse is proportional to the difference 
between the time of occurrence of the threshold crossing and the time 
of the nearest pulse of the displaced sampling wave. The displaced 
sampling wave is 



qd (t) = q(t - T/2) = E 8{t-nT - T/2 - yT), 



(19) 



where y is the phase measured in fractional signal periods. If the axis 
crossing following the with sampling time is displaced by a m T, 



s(mT + T/2 + <x n T) = 0, | a m + 7 I < 4. 



(20) 
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Fig. 1 — Threshold crossing timing recovery system. 
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then the error signal is 

e(t) = E£i(a*-7)5(*-mT), (21) 

m 

where if, is a gain constant. The error pulse has been represented in 
(21) as an ideal impulse function. Since the error signal will be passed 
through a narrow low-pass filter, the response will be virtually identical 
to that when a more realistic pulse of the same area is used. Similarly, 
the effects of variation of the position of the pulse within the interval 
may be neglected, since the low-frequency components of the error 
signal are substantially uneffected. 

We will now determine the threshold crossing time a m as a function 
of the signal overlap and noise. Substituting (1) into (20) yields 

E a m - k f(kT + T/2 - PT + a m T) + n(t) = 0. (22) 

ft 00 

If 

KT/2) + K-T/2) > Z \KkT + T/2)\+n{t) (23) 

then a crossing will occur following the with bit, if and only if a m = —a m+1 . 
When a m = —a m+1 and (23) holds, (22) may be written as 

a m [/(7'/2 - (3T + a m T) - j(-T/2 - 0T + a m T)\ 

+ £ a m . k f(kT + T/2 - $T + a m T) + n(/) = 0. (24) 

**0.-l 

If a m — /3 is small, we may approximate (24) by the first terms of its 
Taylor series expansion. 

aMT/2) - K-T/2)] + o„(a M - P)T[f(T/2) - f(-T/2)] 

+ E a~*KM + T/2) + n(t) « 0. (25) 

Let the abscissa of the function f(t) be adjusted so that 

KT/2) - /(-7V2). (26) 

Then define 

& = 1 f {-T/2) - f(T/2). (27) 

We may now solve (25) for a n . 

a m M (3 + ^ [ E a m . k j{kT + T/2) + n(0]. (28) 
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If we let 

d- fO. if «» = a n+1 (29) 

U, if «n = — a n+ i 
then the error signal (21) becomes 

e{t) = K x t dip - 7 + ^ [ Z a n . k 1(kT + T/2) 

+ n(0]} *(' -«T)- ( 3 °) 

The error signal is passed through the filter (?i(co) and then shifts 
the phase of a reference source. The reference source may be a local 
oscillator whose frequency is tuned as closely as possible to the signaling 
rate. Alternatively, the reference source may be derived from pilot 
tones which are transmitted along with the data. In the latter case, 
there is no error in the average frequency of the reference source, but 
its phase may be poorly related to that of the data signal and may 
also be perturbed by noise. In either case, the reference source generates 
a signal of the form 

r(t) = Z 5(t-nT - rT). (31) 

n = -oo 

When a local oscillator is used, 

t « Aft + 0, AJT « 1 (32) 

where A/ is the frequency offset of the oscillator and 9 is an arbitrary 
constant. When the reference source is derived from pilot tones, 

r = d + 1,(0 (33) 

where y(t) is a zero mean random variable. 

The sampling wave is formed by shifting the phase of the reference 
source by an amount proportional to the value of the filtered error 
signal. 

q(t) = J2 8(t-nT - tT - K 2 vT), (34) 

where v is the filtered version of e(t) and K 2 is the proportionality con- 
stant. Comparing (34) with (2), 

y = t + K 2 v = t+ K 2 J fc(J - x)e(x) dx. (35) 
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In the phase-locked loop type of control system, the frequency of 
the reference source is adjusted rather than its phase. This, however, 
is completely equivalent to integrating the error signal prior to phase 
adjustment. The block diagram is therefore valid for the phase-locked 
loop provided that the filter includes a pole at the origin. As will be 
shown, this pole is highly desirable and, in some cases, absolutely 
essential. 

Substituting (35) into (30) 

e(t) = K x 2 d n \fi - t - K 2 j g,(t - x)e{x) dx 

+ ff [JE a*->W + r/2) + n(0]} S(t - nT). (36) 



Let 



*(0 = ^# giW 



e(0 = 2 e ( n ) *(' _ nT ) 



(37) 
(38) 
(39) 



and normalize the time variable so that T = 1. Then (36) can be 
written as 



e 2 (n) = 2dAfi - r - £ ff 2 (n - k)e 2 (k) 

+ TlZ a n . k 1(k + 1/2) +n(t)]\- (40) 

0*^o.-i ) 

A model of the threshold crossing timing recovery system which con- 
forms with (40) is shown in Fig. 2. This is not a time-invariant linear 
system because of the presence of the multiplier. 
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Fig. 2 — Model of threshold crossing system. 
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In Appendix A, it is shown that 



e 2 (n) = - r - g 2 (0)e 2 (n) - £ g 2 (n - fc)e a (fc). (41) 

*--co 

This system may be readily analyzed either by means of the z-trans- 
form, or equivalently, by discrete Fourier analysis. 2 The discrete 
Fourier transform is given by 

X(a>) = E x(n)exp(-jcon). (42) 

n= — oo 

If G a («) is bandlimited to | « | < r, which is approximately true in 
all cases of interest, then these transforms will coincide with the true 
Fourier transforms. 

The solution of (41) in terms of Fourier transforms is 



rr _ p( u ) - t(«) 

Lz(a,) ~ 1 + fc(0) + QM' 



(43) 



where 



g 2 (0) -j£ftW*». (44) 

The static error can now be found from 

y(u) - /3(co) = G 2 (cS)E 2 (a>) + t(w) - |8(») (45) 

** ~ M ' i + \m + GM *** " ^ (46) 

In particular, if j8(0 is a constant, O , and t(J) is given by (32), then 

to - a. «(«) = 2. 1+ 1 J 'f ) giM B a/ r« + («- M«4] W 
^ - * - i+U)+°U > A " + e ~ M 



do) 



1 + g a (0) 



.1 + ft(0) + G a («). 



(48) 



If G 2 (0) is finite, then the first term is a steadily increasing error 
and the system fails. If A/ = 0, the system does not fail, but a static 
error will be present due to the arbitrary phase values, ft, and 6. It is, 
therefore, highly desirable that G(u>) have a pole at the origin. In the 
presence of frequency offset, this pole is essential. In this case, only 
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the last term in (48) is nonzero, so that a frequency offset leads to a 
static error. The system is completely insensitive to the values of the 
arbitrary phases, /3 and 0, except during initial start-up. 

For the system to be stable, it is required that 1 + g 2 (0) + (7 2 (co) 
have no zeros in the half plane Im(aj) < 0. In most cases, g 2 (0) « 1, 
so that the usual criteria for stability apply to G 2 (<a). 

We now wish to calculate the variance of 7, or the mean square 
jitter. From (40) and (41), and assuming the static error to be constant, 

e 3 (n) = e 2 (n) — e 2 = 2d„\ e 2 + g 2 (0)e 2 - Z Q*X n ~ fye 3 (k) 

- v(t)+jn(t)j+z(n) -7 2 , (49) 



where 



Let 



z(n) = \d n a n £ a„_J(fc + 1/2). (50) 

hno.-l 



x(n) = 2d]^e 2 + g 2 (0)c 2 + jn(t)j + z(n) - e 2 . (51) 

Then (49) may be written as 

e 3 (n) = x(n) - 2d nV (0 ~ 2d n £ g 2 (n - k)e 3 (k). (52) 

The zero mean component of the output phase error is 

7>fa) - y(n) - 7 = 77(0 + Z 9^n ~ *)*(*) (53) 

7i(w) = 77(0 + Z g 2 (n - k)[x(k) - 2d k7l (k)]. (54) 

The autocorrelation of (54) is 

Z Y,E{[2d k g{n - k) + 8 nk )[2d l9 (n + m - I) + 3 n+m .<]7*7<} 

= Rr,(m) + Z Z 0»(» - k)g 2 (n + m - l)R z (k - T), (55) 

k I 

where E(v) denotes the expected value of v and R v (m) = E[v(n)v(n + m)] 
is the autocorrelation of v. 
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In almost all cases of practical interest, g 2 (0) « 1. Then we may 
approximate 

* = 2 (O)fc & (56) 



d n y(ri) = \z « (57) 

and 

* = /?~L « ft - 7. (58) 

1 + 02(0) 

Subject to these approximations, we can evaluate the discrete Fourier 
transform of (55). After some algebraic manipulation, the result ob- 
tained is 

| 1 + G 2 (o>) | 2 S y (o) + | G 2 (u) | 2 R y (0) = £» + f G,(«) | 2 &(o>), (59) 
where 

CO 

£.(«) = £ #.("0 exp (-;wco) (60) 

m = — oo 

is the power spectral density of v. 
The variance of y is calculated as 

<r\ = R y (0) = - f S y (a>) do) (61) 

7T Jo 

i r I g 2 (u) |» i r g^) 

2 = 7T Jp | 1 + O2M I X Jp I 1 + Cr 2 (fa)) I 

■K J I 1 + Cr 2 (w) I 

Since (? 2 (w) is narrowband, we may assume that 



?2 aw 



(62) 



GM 



+ <?*(«) 



dw « 1 (63) 



r - 

io I 1 
and therefore, 

^^tJo I 1 + G 2 (o) «»*•+;./, |i +(?,(«) |- *"• (64) 

The second term indicates that low-frequency components of the 
reference source noise are attenuated while high-frequency components 
are not. Therefore, if the reference source is derived from transmitted 
pilot tones, it should be filtered to a narrow bandwidth before being 
used in the timing recovery system. 
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In Appendix B, &(&>) is evaluated for | w | « 1. This is the only region 
of interest when G 2 (o)) is a sufficiently narrow low-pass filter. 

&(«) = e\ +p [Rn(0) - R n (l)] + A 1 + A 2 co 2 , (65) 

where 

4:4i; w + 1/2) - /(* - 1/2)] 

O k CO 

• [2/(fc + 1/2) +/(-*+ 1/2) - f(-k - 1/2)] (66) 
and 

A > = 2P {~ 4 ^ 1/2 ) + 2 4 ?„ /( * + 1/2)/(fc " 1/2) 

- Z fc 2 [/(^ + l/2)-/(fc-l/2)][/(-A;-M/2)-/(-fc-l/2)}- (67) 



k--oo 

If we let 



and 



i r I 

Jl = ~ l 7 



G 2 («) 



CJ 2 



1 r » 

= - / w T 

7T Jo I - 1 - 



+ <7 2 (a>) 

G 2 (a>) 



dco (68) 



+ G 2 (<o) 



d<o (69) 



then (64) becomes 



a* = ? 2(ai + p [8.(0) - «„(!)]«, + A,», + A 2 col 



+ i P I i If) x ■» *•• (70) 

f JO 1 + Cr 2 (w) | 



This equation is given in the summary with the normalization T = 1 
removed. An application to a typical data transmission system is given 
in Section IV. 

It should be noted from (69) that a> 2 will be unbounded unless G s {<a) 
has at least two more poles than zeros. Good design of G 2 (u) requires 
that the second pole (assuming no zeros) occur somewhere in the 
vicinity of gain crossover. In this case, o> 2 is approximately equal to w, . 

The first term of (70) indicates that the standard deviation of the 
jitter caused by frequency offset will be much less than the mean 
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of that error. The second term is the jitter produced by the additive 
noise. In typical data transmission systems, the signal is filtered such 
that R„(l) ~ 0. In that case, the variance of the jitter is directly pro- 
portional to the noise power. Of particular interest is the jitter produced 
by signal overlap. It can be seen from (66) that A x = if j(t) is an even 
function. In that case, the jitter variance is proportional to the cube 
of the system bandwidth, and can be made quite small by the use 
of narrow filtering. If f(t) is not an even function, then the third term 
will usually be much larger than the fourth term. In either case, the 
jitter will greatly increase as the filter used to shape f(t) is made sharper. 

III. THE SAMPLED-DERIVATIVE DETECTOR 

An alternative method of adjusting the phase of the sampling wave 
makes use of the time derivative of the signal at the sampling times. 
Implementation of a sampled-derivative timing recovery system is 
about equally complex as a threshold crossing system. 

Except for the manner in which the error signals are generated, 
the control loop is the same for both systems. Fig. 1 may be used to 
describe the sampled-derivative system if the delay in the feedback 
path is eliminated and a sampled-derivative detector is substituted 
for the threshold crossing detector. 

The sampled-derivative detector generates an error pulse during each 
bit interval whose amplitude is proportional to the time derivative 
of the data signal at the sampling time, multiplied by the polarity of 
the signal at that time 

e(t) = K 3 2 sgn [s(?iT + y T)]s'(nT + yT) S(t - nT). (71) 

However, the output data is generated by setting 

& n = sgn [s(nT + yT)], (72) 

where d n is the receiver decision on the nth bit. If the error rate is low, 
&n = o-n with high probability, and (71) may be approximated by 

e{t) ^ # 3 Z a n s'(nT + yT) 8(t - nT) (73) 

if the effect of errors is neglected. 
Using (1), 

e{t) = K 3 J2 aJE a n -. k f(kT + yT - f3T) + n'(t)} 8(t - nT). (74) 
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The abscissa of j{t) in this case is adjusted such that the origin coin- 
cides with the peak of j(t). 

/'(0) = 0. (75) 

If the phase error 7 — /3 is small, we may approximate (74) by the 
first terms of its Taylor series expansion 

e(f)ttK* E a n [a n ( 7 -0)77"(O)+ E <*»-*/ W +*'(')] 6(t-nT). (76) 

As in the previous case, 

e(t) = E e(n) 8(t - nT) (77) 

n 

y = K 2 f 9l (t - x)e(x) dx + r (78) 

and we normalize the time variable by setting T = 1. Equation (76) 
may now be written as 

e(n) = K 3 [(t - /3)/"(0) + a n E o._ 4 /'(fc) + a»n'(0 



+ X2HO) E gi (n-k)e(k) 1. 



(79) 



Let 



«»(0 - "IF^T (80) 



and 

*(*) = -K 2 K 3 1"(0) gi (t). (81) 

Then (79) becomes 

e 3 (n) - |8 - r - y(n) - E £ 3 (n - fc)c 3 (fc), (82) 

where 

»(») = TTfer [E «--*/'(*) + »'»]. (83) 

Unlike the threshold crossing system, the sampled-derivative system 
is a time-invariant linear one when the phase error is small. A model 
of the system conforming with (82) is shown in Fig. 3. This model 
may be readily analyzed because of its linearity. 
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fi-y 



O 



G 3 (o>) 



Fig. 3 — Model of sampled-derivative system. 

n 

y(n) = 2 Oafa - k)e 3 (k) + r 

7(n) = E g 3 (n - k)[fi - y(k) - y(k)] + r. 

t--00 

The mean error may be found from 

7(0) [1 + G 3 (a)] = <?,(«)jS(«) + r(u>) 

t(co) - /3(o>) 



7 (o>) - /3(co) = -£,(«) = 



1 + G 3 (w) 



(84) 
(85) 

(86) 
(87) 



Equation (87) is identical to (46) if 62(a) is substituted for(7 2 (o>)/[l + 
<7 2 (0)]. All the comments of Section II concerning the static error and 
the desirability of G(u) having a pole at the origin therefore, also apply 
to the sampled-derivative system. 

The variance of 7 when the mean error is constant will now be found. 

7.(») =y(n) ~y = -Z g 3 (n - k)[y(k) + 7,(*)] + n(0, (88) 

k--n 

where ij (I) is again the reference source jitter. In terms of the power 
spectral densities, 



s» = 



G 3 (o)) 



1 + G-M 



8M + n + GMr 



In Appendix C it is shown that, for | w | « 1, 
where 



(89) 



(90) 



A * = P(0) .§. f'(*W<*> + ^"*>1 



(91) 
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and 



"4 — ctx'f 



2/"'(0) ut 



Z k 2 f(km-k). 



The variance of 7 is then 



r y = — J S y (co) dw 

7T Jo 

fi^(o) , , . , 3 , 1 r &,(«) , 

7 (0) 7T J I 1 + 0- 3 (co) J 



where 



and 



1 r 

T Jo I 1 

3 1 r 2 

w 2 = - I w 

7T Jo 



O a (fa) 



+ G 3 (o;) 



G 3 (co) 
1 + C7,(«) 



c/w. 



(92) 

(93) 
(94) 

(95) 

(96) 



There is a very strong similarity between (94) and (70). The last 
terms are identical, so that it is just as important in the sampled- 
derivative system as in the threshold crossing system that high-fre- 
quency noise components be removed from the referenced source prior 
to use for timing recovery. 

The jitter due to additive noise is proportional to the power of the 
derivative of the noise. The example in the next section illustrates 
that this is not serious if the noise is bandlimited to the same frequency 
range as the signal. However, if any high-frequency noise is allowed to 
enter the receiver beyond the signal filter, the jitter will be greatly 
increased. 

The jitter due to signal overlap is very similar to that of the threshold 
crossing system. If f(t) is an even function, then its derivative will 
be odd, and A 3 = 0. Both 4 3 and A t increase markedly as the spectrum 
of /(/) is made sharper. 

Finally, unlike the threshold crossing system, there is no additional 
jitter term due to static phase error. This jitter component is eliminated 
because there is an error pulse generated during each bit interval. 

IV. AN EXAMPLE 

In order to illustrate the results of the previous two sections, the 
output jitter of both a threshold crossing timing recovery system and 
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a sampled-derivative system will be calculated for a typical received 
data transmission signal. 

The signal to be considered is one using half raised-cosine amplitude 
shaping which is distorted by linear delay distortion. For simplicity, 
it is normalized so that T = 1 and the undistorted peak of the signal is 1. 



F(u) = A(cS) exp [#>(«)] 



(97) 



(98) 



(99) 



The outline of the "eye pattern" for this signal is shown in Fig. 4, 
along with the central portion of f(t). The eye pattern is formed by 
superimposing the signals of all possible message sequences. Closing 
of the eye is due to signal overlap. 

The time of maximum eye opening is the optimum sampling time in a 
minimax sense. When such a sampling instant is chosen, then the error 





1, 


»<«<! 


A(a) = ■ 


2 W ~2 
cos 2 , 


IT . ."&TC 

2 < w < T 


*>(<*) = 4^<" 2 . 


u > 0. 




~n r 

SD TC IM 

Fig. 4 — Eye pattern outline and mean sampling times for a typical distorted 
data signal. 
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probability in the presence of additive noise is minimized for the most 
adverse message sequence, where all adjacent bits overlap in a manner 
that subtracts from the bit being detected. 

The mean sampling times are also shown in Fig. 4 for the case of 
no static error. The mean sampling time of threshold crossing system, 
TC, is such that j(TC + J) - J(TC - |). The mean sampling time 
of the sampled-derivative system, SD, coincides with the peak of jit). 
It is seen that the threshold crossing system chooses a better mean 
sampling time than the sampled-derivative system, yet both systems 
miss the maximum of the eye opening by a large amount. 

In order to compare the jitter due to noise, a particular noise spec- 
trum must be considered. Here it will be assumed that the noise is 
white noise which has been passed through a receive filter matched 
to the undistorted signal. In this case, the noise power spectral density 
is of the form 

m - 1 £"* (ioo) 

- / Ado) dco 

IT Jo 

so that 

R n (0) = vl , (101) 

RnO) = 0, (102) 

and 

/ u A(co) dco 

R'n'(0) = —^ cr; . (103) 

/ A (to) dco 

Jo 

From (70), the rms jitter due to noise for the threshold crossing 
system is calculated to be 



<r y = 0.615o- B V«i • (104) 

For the sampled-derivative system, this quantity is calculated from (94). 

a y = 0.612a„V^ . (105) 

The results are virtually identical. In either case, if a n = 0.1 (signal- 
to-noise ratio of 20 dB) and w, = 0.01, then the rms jitter due to noise 
alone will be 0.01 percent. It should be mentioned that several other 
signal pulse shapes were examined, and it was found that the jitter 
due to noise was not very sensitive to pulse shape. 
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In order to calculate the jitter due to signal overlap, A x and A 3 
were computed from (66) and (91). A 2 and A 4 were small enough to 
have negligible effect on the jitter. The resultant rms jitter is 0.287 Vw] 
for the threshold crossing system and 0.265 \Zcol for the sampled- 
derivative system. If w, = 0.01, the jitter is 2.87 percent and 2.65 
percent, respectively. This is by no means negligible, and illustrates 
the need for very narrow filtering in the timing recovery control loop. 
Again there is little difference in the performance of the two systems. 

To observe the effects of asymmetry of the signal pulse, let us con- 
sider the same signal without phase distortion. Both timing recovery 
systems will then set a mean sampling time at the best point. The 
jitter due to signal overlap is greatly reduced since A x and A 3 are zero. 
The computed values of A 2 and A A are 0.11 and 0.32, respectively. 
If o) 3 = 0.01, then the rms jitter due to signal overlap is only 0.01 
percent for the threshold crossing detector and 0.03 percent for the 
sampled-derivative detector. Both values are completely negligible. It 
may be concluded from this calculation that both timing recovery 
systems are very sensitive to asymmetry of the signal waveform, both 
in terms of choosing the average sampling time and the resultant jitter 
about that time. 

V. AN IMPROVED TIMING RECOVERY SYSTEM 

It was seen in the previous example that both the threshold crossing 
timing recovery system and the sampled-derivative system led to 
average sampling times which differed considerably from the time of 
maximum eye opening. However, it is possible to modify the sampled- 
derivative system so that it does seek the time of maximum eye opening 
as the average sampling time. 

At any time t , the signal amplitude for the worst message sequence, 
assuming the current bit is 1, is 

D(t ) = KQ - E | 1(t + kT)\. (106) 

In the region where the eye is open, D(t ) > 0, and the eye opening 
is equal to 2D(t ). If a sampling time t is chosen such that D(t ) < 0, 
then errors will occur for some sequences even in the absence of noise. 

An experimental examination of the eye patterns of a large number 
of actual data transmission systems indicates that D(t ) is almost 
always a concave function of / • Therefore, if t is adjusted according 
to the gradient of D(t ), then the maximum of D will be found. 

It is therefore desired to generate an error signal whose average 
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value is 

e(5 = KD'(t - 0T + yT). 
If we again normalize the system so that T = 1, 



(107) 



e{t) = K[V{t - + 7) 

- £ 1V + * - + 7) sgn /(/ + k - + 7)]. (108) 

Equation (108) exists and is continuous except at those points t k 
where j(t k + k - + y) = 0. 

Fig. 5 is a block diagram of a system which generates an error signal 
whose average is given by (108). The first term of (108) is the average 
error signal of the sampled-derivative detector discussed in Section III. 
The improved timing recovery system therefore, will consist of a 
sampled-derivative system with added correction signals. Enough cor- 
rection terms are used to account for those adjacent bits which may 
be expected to overlap significantly into the bit interval under con- 
sideration. 



SHIFT REGISTER 




Fig. 5 — Improved timing recovery system. 
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The kth correction signal must have an average value 

7M = K1'(k + sgn f(k + ft), (109) 

where 

h = i - P + 7- (HO) 

This correction signal is formed by first generating an auxiliary 
signal x k (ri) whose polarity is expected to agree with that of f(k -f- J x ). 
The derivative of the signal at the current sampling time is multiplied 
by the polarity of the signal at a time displaced from the current time 
by k bit intervals. The result is multiplied by the polarity of the auxiliary 
signal to form the correction signal. 

e k (n + U) = Ks'(n + h) sgn s(n — k -f /,) sgn x k (n + *,). (Ill) 

In order to account for overlap into leading pulses as well as lagging 
pulses, a fixed delay must be built into the system, as indicated in 
Fig. 5. This delay is equal to half the shift register length, so that the 
central cell of the shift register stores the polarity of the current bit, 
while the other cells store the polarities of preceding and succeeding bits. 

The auxiliary signal is formed by multiplying the value of the signal 
at the current sampling time by the polarity of the signal which pre- 
ceded this signal by k bit intervals. If k is negative, the polarity of a 
succeeding bit is used. The resultant is filtered by a narrow filter, 
H(o)), to form the auxiliary signal x k . 

x *( n ) = Z^( n — w)s(m) sgn s(w? — k), (112) 

where the time displacement t t is ignored. 

If we assume that the error rate is low, as was done in Section III, 
then we may approximate sgn s(n) = a n with little loss of accuracy. 
Using this approximation and substituting (1) into (112), 

x k (n) = Z h(n - m)[a m - k n(t) + f(k) + a m _* £ a m .J(j>)]. (113) 

The mean of x h is 

x~ k = j(k)H(0). (114) 

In Appendix D it is shown that the variance of x k is 

<£ - ;f f I //(«) I 2 °l + f P I F(v) T *• + Pub) <*«. (H5) 
where P k (oi) is the Fourier transform of 

P*(0 = /(* - 0/(* + 0- (116) 
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If H(co) is sufficiently narrowband, then it can be seen from (115) 
that the variance of x k will be small. Then, if f(k) is not too small, 
it may be expected that the polarity of x k will agree with that of j(k) 
with high probability. 

We will now examine the mean value of the correction signal. Using 
the assumption of low error probability and substituting (1) into (111), 
the correction signal is 

e k (n + /,) = /Co„_, £ o.-i/'(» + k) sg" x k (n + /,). (117) 

i 

In order to rind the mean of (117), we make the approximation that 
x k (n) is independent of a n _* and a n _, . As a justification of this ap- 
proximation, note from (113) that 

^ | a„_* , a„_, = Wc)H(0) + f(l)h(0) + f(2k - l)h(l - k) (118) 

x k | a„_ t , a„_, « .r* (119) 

since /i(?i) <<C #(0) for a narrowband filter. Let 

P 4 = Prob [sgn x k = sgn f(k)]. (120) 

Then 

el « K(2P - l)j'(k + /,) sgn /(/c + /,)• (121) 

When the magnitude of /(/c + /,) Ls sufficiently large, P tt 1 and the 
mean of the /cth correction signal is approximately the desired value. 
In the vicinity of a zero of f(k + t t ), | < P < 1, and the correction 
signal will at least have the correct polarity, although not the correct 
magnitude. Ideally, the correction term should be a discontinuous 
function of t x at a zero of f(k -f- ti). The actual correction signal will 
have a mean value which is continuous, but the sharpness of change in 
the vicinity of a zero will increase as H(o)) is made narrower. 

The rms jitter of this timing recovery system is extremely difficult 
to evaluate because of the presence of many nonlinear operations. 
However, this jitter may be expected to be much greater than that 
of a sampled-derivative system, since each correction signal may be 
expected to introduce jitter of the same order of magnitude as the 
main error signal. Narrow filtering in the control loop is therefore 
essential. 

The mean sampling for the example of Section IV when this improved 
timing recovery system is used is shown in Fig. 4 as "IM". It is seen 
that the time of maximum eye opening has been found. For this example, 
only one leading and one lagging correction term were sufficient to 
choose this mean sampling time. 
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VI. EXTENSIONS 

Each of the timing recovery systems described here may be modified 
to work with ??i-level digital data signals instead of only binary formats. 
Since in multilevel systems the eye is much narrower, the effects of 
static timing error and jitter are much more serious. 

A threshold crossing detector may be constructed which generates 
an error pulse whenever the signal crosses any one of the m — 1 thresh- 
olds. During any signal interval, any number of error pulses between 
zero and m — 1 may be present. Such a system is extremely difficult 
to analyze, but has been found to work well in practice, provided 
that some auxiliary means is used to correct the mean sampling time. 

The sampled-derivative detector is very easily extended to multi- 
level systems if the signal derivative is multiplied by the output symbol 
value in forming the error signal. If the a„'s are scaled so as to form 
a set of unit variance, then the analysis of Section III applies directly. 

The improved timing recovery system is modified for use with multi- 
level signaling in a manner similar to that of the sampled-derivative 
system. However, the signal margin against noise for the worst message 
sequence is now 

D(Q = KQ - (m - 1) Z I K*o + kT) |. (122) 

Comparing this criterion with that of (106), it is seen that each of the 
correction signals must be weighted by the quantity m — 1 in order 
to find the time of maximum eye opening. 

Extension of these techniques to partial response systems is also 
straightforward. Since the modifications depend on the particular 
partial response system used, a description will not be presented here. 

All of the systems analyzed here used linear control loops. This 
permitted the calculation of jitter variance in terms of loop bandwidth. 
However, the implementation of these systems may frequently be 
simplified considerably by using nonlinear control systems. A particular 
method which has met with practical success uses the polarity of each 
error pulse to adjust the sampling phase by a small fixed increment. 

APPENDIX A 

Evaluation of e 2 * 

Equation (40) is of the form 

e 2 (n) = 2d\c n - Z 9*(n ~ k)e 2 (k) , (123) 



* The approach used here was suggested by J. E. Mazo. 
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where d n = or 1 with equal probability and are independent, The 
random variables c n are uncorrected with the d„'s. 
From (123) and the properties of d„ , 

e 2 (n) = d n e 2 (n) (124) 

and 



d n e 2 (k) = \e 2 (k), if n > k. (125) 

We first find the average of (123) over all c„ and all d k , k < n. This 
average of a random variable v will be denoted by (v), while the average 
over all c„ and d n will be shown as v. 

(e 2 (n)) = 2d n 7 n -2E g 2 (n - k)(d ri e 2 (k)). (126) 

Using (124) 

(e. 2 (n)) = 2rf,A - 2g 2 (0)(e 2 (n)) - 2 E g 2 (n - k)(d n e 2 (k)). (127) 



The only random variable in (127) is d„ . The overall average, e 2 {n), 
is therefore the average of (127) over d„ . Using (125), we obtain 



e 2 (n) = c„ - 2g 2 (0)e 2 (n) - £ g 2 (n - k)e 2 (k) (128) 



e 2 (n) = c n - g 2 (0)e.M - £ oM ~ *)«•(*) ( 129 ) 

it— BO 

which is the result shown in (41). 

APPENDIX B 

Evaluation oj S z (w) 

From the definition of d„ in (29), we may express 2d n as 

2d n = 1 - ba„. (130) 

Then (51) may be rewritten as 

a:(n) w -m.*A + «(«) + | («, - «-..,)n(/), (131) 

where the term g 2 (0)e 2 has been neglected. 

We wish to evaluate the power spectral density of x. It will first be 
shown that the approximation of x given in (131) is zero-mean. Since 
the a n 's are zero-mean and independent, and the noise n(t) is zero- 
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mean and independent of all other random processes, then the first 
and last terms of (131) are zero-mean. The mean of z is readily found 
after substituting (130) into (50). 

z(n) = I (a. - a n+1 ) £ «„_./(/,- + J) (132) 

zfyt) = 0. (133) 

The mean value of x(n) as given in (131) is therefore zero. It may 
similarly be shown that the three terms of (131) are mutually uncor- 
rected. The autocorrelation of x is then 



R z (m) = a n a n+1 a n+m a n+m+ i ei + i£ 2 (ra) 

1 



+ p (a n — a n+1 )(a n+m — a n+m+ ,)R n (m) (134) 

RM = 3 + R,(0) + pfi»(0) (135) 

R x (±l) = fl,(±l) - Afl„(±l) (136) 

jR.(m) = fl.(m), m ^ 0, p* 1. (137) 
From (132), 

Rz{m) = T2 2Z 5Z (°» ~ a»+i)(o»+« - On+m+i)a»-*a»+—i 

**(>,- 1 J 4=0.-1 

/(ft + fi/(l + J) (138) 

RM = 1 Z ft* + h) (139) 

O /M0.-1 

fl.(±l) - -A [ Z /(* + §)/(* - I) + /(-!)/(*)] (HO) 

R,(±m) - £ [/(m + |) - ft* - DM-« + i) - K-m - J)], 

w ^ 0, ^ 1. (141) 

The power spectral density of x can now be calculated for \ ta \ < r. 

SM = J2 R*(m) exp (-jum) (142) 

m = — » 

CO 

S x (u) = R x (0) + 2 £ik(m) coswm. (143) 
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Since x will be passed through a very narrow filter, we are interested 
in the spectrum of x only in the region | o> | <K 1. In this region, (143) 
may be approximated by 

CO 00 

&(«) « R x (0) + 2 £ fl x (m) - a. 2 £ m 2 R s (m). (144) 

m=l m = l 

This approximation is valid if the third derivative of S x (u) is bounded. 6 
This will be true if R z (m) decreases as 0(m~ 5 ). If the Fourier transform 
of j{t) is continuous, then /(m) decreases as 0(m~ 2 ). In this case, it 
can be seen from (137) and (141) that R x (m) decreases as 0(m~°), 
so that the approximation (144) is valid for | u | « 1. 

S,(«) « ? a + |b |/2.(0) - fi„(l) + Z f (& + §) 

- Z fQo + *)/(* -h) - /(-§)/(!) 

«• * ii.ii 

+ E [/(m + *) - Km - MK-m + *) - /(-m - J)]} 

rn-2 ' 

+ « 2 { 2 w + *)/(* - i) + /(-§)/(§) 

- Z m a [/(m + ft) - ftm - ft][/(-m - *) - /(-m - *)]}■ (145) 
After some manipulation, and using (26), (145) may be reduced to 
SM = el + 1 |#„(0) - B„(l) + | jL [/(* + I) ~ /(* " *>] 

• [2/(/c + *) + /(-* + »-/(-*- *)]} 
+ ^ {-4f (J) + 2 JC f(A + |)/(* - $ 

- E fc 2 [/(* + |) - f(k - h)M-k + ft - /(-* - *)]V (146) 



APPENDIX C 



Evaluation oj *S„(a>) 

We wish to find the power spectral density of 



/ (y) k#0 
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The autocorrelation of y is 

R*(in) = 77727^- [£ £ a„a n+m a„- k a n +„-,1'(k)j'(£) 
! \\>) k*n i*u 



+ a n a n > m n'(t)n(t + m)\ (148) 

RM = tttLt [Z f*<*) - «'<©]■ (149) 

7 w! **o 



For m t± 0, 



^ (m) = rko) r(m)r( " w) ■ (150) 



Under the same conditions stated in Appendix B, the power spectral 
density of y may be approximated in the region | u | <3C 1 by 

S y (a>) « fi,(0) + 2£ #„(w) - co-' £ m 2 R v (m). (151) 

Using (75) we obtain 

SM M p^y {-^'(0) + JC f <*)IF(*) + f(-*)] 

- K £ *V'(*)f(-*)[- (152) 

APPENDIX D 

Evaluation of a\ k 

The variance of .r* is the mean square value of the zero-mean random 
process 

x k (n) - x = Y,h(n- m)a m - k [n(t) + £ a»-j(p)l (163) 

Since the noise and the message are independent, the variances 
of the two components of (153) will add to form the total variance. 
The variance of x k due to the noise is 



°i = £ £ Hn ~ in)h(n — q)a n - k a Q - k n{t m )n(t Q ) (154) 

m q 

** - £ ^fa ~ ™> 2 (/) , 155) 

„* = -L a l [ I #( w ) I 2 rf w< (15(}) 
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The total variance of x k is then 



<rli- = it; + E E E Z h(n - >n)h(n - q)a m - k a v - k a n - p a Q - r 

m >i p*k r * k 

■MK'i). (157) 

Nonzero contributions arise from those terms where q = m and 
r = p and from those terms where r = 2k — p and q = m + k — p. 

vlu = ffl + E E [h*(n - m)f(p) 

+ h(n - m)h(n - m - k + p)f(p)f(2k - p)\. (158) 

Any one of the terms of (158) is small compared to the sum. We 
may therefore approximate (158) by including the missing p = k terms. 

trl k ft* a\ + A I | H(u) \ 2 du J | F(cS) | 2 dta + d , (159) 

where 

°l = E Z A(m)*(i« - * + P)1(p)1(2k - p) (160) 

m p 

and it is assumed that H(u>) is bandlimited to | o> | < x. 
a 2 = ~ JJff H{o>)H{u)F(v)F(y) £ exp [jm(co + u)] 

• ^ exp [jp(w + r — //)] exp [jk(2y — u)] du du dv dy (161) 

V 

*\ = ^ ffff H(u)H(u)F(v)F(y) 5(co - u) 6(u - v - y) 

•exp [jk(2y - u)] du> du dv dy (162) 
<r\ = A r H(f*)H{-<*) r/co / F(y - »)f(y) exp [j7,-(2v/ - «)] (fy (163) 

o-r: = A f* I //(«) | 2 rfw /" F*(" - //) exp [-#(« - y)] 

■F(y)exp(jky)dy. (164) 

The second integral in (164) may be recognized as 2t times the 
Fourier transform of 

p k (t) = j(k - /)/(/,- + /) (165) 
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so that 

al = ± T \ H(f*) | 2 P*(«) da. (166) 

ZT J-oo 

Substituting (166) and (156) into (159) 

* = h r ' H((S) |2 [*■ + hfj F(u) |2 du + Pi(w) ] *•■ (167) 
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